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ABSTRACT 

HH 212 is a highly coUimated jet discovered in H2 powered by a young Class source, IRAS 05413- 
0104, in the L1630 cloud of Orion. We have mapped around it in 1.33 mm continuum, ^^CO ( J = 2 — 1), 
13C0 (J = 2 - 1), C^^O (J = 2 - 1), and SO {Jk = 65 - 64) emission at - 2.5" resolution with the 
Submillimeter Array. A dust core is seen in the continuum around the source. A flattened envelope is 
seen in C^^O around the source in the equator perpendicular to the jet axis, with its inner part seen in 
^3 CO. The structure and kinematics of the envelope can be roughly reproduced by a simple edge-on disk 
model with both infall and rotation. In this model, the density of the disk is assumed to have a power-law 
index of p = —1.5 or —2, as found in other low-mass envelopes. The envelope seems dynamically infalling 
toward the source with slow rotation because the kinematics is found to be roughly consistent with a 
free fall toward the source plus a rotation of a constant specific angular momentum. A ^^CO outflow is 
seen surrounding the H2 jet, with a narrow waist around the source. Jetlike structures are also seen in 
"'^^CO near the source aligned with the H2 jet at high velocities. The morphological relationship between 
the H2 jet and the ^^CO outflow, and the kinematics of the ^^CO outflow along the jet axis are both 
consistent with those seen in a jet-driven bow shock model. SO emission is seen around the source and 
the H2 knotty shocks in the south, tracing shocked emission around them. 

Subject headings: stars: formation — ISM: individual: HH 212 — ISM: jets and outflows. 



1. INTRODUCTION 

Despite recent progress, the detailed physical processes 
in the early stages of low-mass star formation are still un- 
certain. Currently, the most detailed and successful model 
of low-mass star formation is the four-stage model pro- 
posed in the late 80's by Shu, Adams, & Lizano (1987). 
In the flrst stage, slowly rotating molecular cloud cores 
form within a molecular cloud as magnetic and turbulent 
support is lost through ambipolar diffusion. In the second 
stage, a molecular cloud core becomes dynamically unsta- 
ble and collapses from the inside-out, forming a protostar 
surrounded by an accretion disk deeply embedded within 
an infalling envelope of dust and gas. In the third stage, 
as the protostar continues accreting mass through the ac- 
cretion disk, a stellar wind breaks out along the rotational 
axis of the system, creating a bipolar outflow. In the fourth 
stage, the infall and accretion terminate, leaving a newly 
formed star surrounded by a circumstellar disk. In this 
paper, we investigate the protostellar collapse and outflow 
in the second and third stages, with the Submillimeter Ar- 
ray (SMA)^ observations around a young protostar, IRAS 
05413-0104, and its remarkable jet HH 212. 

IRAS 05413-0104 is a cold, low-luminosity (~ 14 L©), 
low-mass source (Zinnecker et al. 1992), located at a dis- 
tance of 460 pc in the L1630 cloud of Orion. It is sur- 
rounded by a cold (~ 14 K), flattened (with an aspect ratio 



of 2:1) rotating NH3 envelope with a diameter of 12,000 
AU (Wiseman et al. 2001). It is detected also at 1.1 mm 
(Zinnecker et al. 1992) and 1.3 mm (Chini et al. 1997), 
with a high ratio of millimeter-wave luminosity to bolo- 
metric luminosity. It is thus classifled as a Class source, 
which has much material to be accreted still from the sur- 
rounding molecular gas envelope (Andre et al. 2000). 

HH 212 is a remarkable jet with a length of ~ 240 arcsecs 
(~ 0.6 pc) powered by the IRAS source. It was discovered 
in the rotational- vibrational 1; = 1-0 S(l) emission line of 
H2 at A = 2.122 /zm (Zinnecker, McCaughrean, & Rayner 
1998). It is highly coUimated and symmetric with matched 
pairs of shock knots (knotty shocks) and bow shocks on 
either side of the IRAS source. Strong water masers are 
seen moving at 60 km s~^ along the jet axis near the 
IRAS source (Claussen et al. 1998). Shocked SiO emis- 
sion is seen along the jet axis (Chapman et al. 2002) and 
around the source (Gibb et al. 2004). CO outflow is also 
seen surrounding the H2 jet (Lee et al. 2000). Deep obser- 
vations in H2 with the ESO Very Large Telescope (VLT) 
show a pair of diffuse nebulae near the bases of the jet, 
probably tracing the outflow cavity walls illuminated by 
the bright knotty shocks around the source (McCaugh- 
rean et al. 2002). Based on the relative magnitude of the 
proper motions and radial velocities of the water masers, 
the CO outflow is believed to lie within 5° of the plane 
of the sky (Claussen et al. 1998), and is thus an excellent 
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system to investigate the transverse kinematics of outflow. 

Here, we study in detail the circumstellar envelope 
around the IRAS source, the molecular outflow surround- 
ing the H2 jet, and the interaction between them, with 
the SMA observations in 1.33 mm continuum, ^^CO (J = 
2 - 1), (J = 2 - 1), Ci^O (J = 2 - 1), and SO 

{Jk = 65 — 54) emission at high angular and velocity res- 
olutions. It is currently believed that 1.33 mm contin- 
uum mainly traces dust core, ^^CO (particularly in inter- 
ferometric; observations) mainly traces molecular outflow, 
C^^O mainly traces envelope, ^'^CO traces both molecular 
outflow and envelope, and SO traces mainly shock inter- 
actions. 

2. OBSERVATIONS 

Observations in 1.33 mm continuum, ^^CO (J = 2 - 1), 
13C0 (J = 2-1), C^^O (J = 2-1), and SO {Jk = 65-64) 
emission toward the IRAS source and the H2 jet were made 
between 2004 November and 2005 March on top of Mauna 
Kea with the SMA in the compact configuration (see Table 
1). Seven antennas were used in the array, giving baselines 
with projected lengths ranging from 14.1 to 136 m, result- 
ing in a synthesized beam (with natural weighting) with a 
size of 2.8" X 2.4" at 230 GHz. With a shortest projected 
baseline of 14.1 m, our observations were insensitive to 
structures more extended than 15" (^ 7000 AU) at the 
10% level. The FWHP of the primary beam is 55", and 9 
pointings were used to map the whole jet system. The dig- 
ital correlator was set up with a bandwidth of 104 MHz for 
each band (or chunk). We used 512 spectral channels for 
i^CO, 0^*0, and ^^CO and 256 spectral channels for SO, 
resulting in a velocity resolution per channel of ~ 0.27 km 
s-i for 12C0, C^^O, and ^^CO, and 0.55 km s"! for SO. 
The 1.33 mm (or 225 GHz) contimmm emission was also 
recorded with a total bandwidth of ~ 4 GHz. In our obser- 
vation, a phase cahbrator (Quasar 0423-013) was observed 
every 20 minutes to calibrate the phases of the source. A 
planet (Uranus) was observed to calibrate the fiuxes. The 
data were calibrated with the MIR package adapted for 
the SMA. The cahbrated data were processed with the 
MIRIAD package. The dirty maps that were produced 
from the calibrated data and corrected for the primary 
beam attenuation were deconvolved using the Steer clean 
method. The final channel maps were obtained by con- 
volving the deconvolved maps with a synthesized (Gaus- 
sian) beam fitted to the main lobe of the dirty beam. The 
velocities of the channel maps are LSR. 

3. RESULTS 

In the following, the IRAS source is assumed to have 
a refined position at q:(2ooo) = 05^^43^51^408, (5(2ooo) = 
-01°02'53.13", which was found at A = 3.5 cm with the 
VLA at an angular resolution of 0.3" (Galvan-Madrid et 
al. 2004). The systemic velocity in this region was as- 
sumed to be 1.8 km s"-'^ in CO J=l-0 observations (Lee et 
al. 2000) and 1.6 km s~^ in NH3 observations (Wiseman 
et al. 2001). Here, it is assumed to be 1.7 ±0.1 km s~^. 
Throughout this paper, the observed velocity is the veloc- 
ity with respect to this systemic velocity. For comparison, 
our observations are always plotted with the H2 image of 
the jet recently obtained with the VLT (McCaughrean et 
al. 2002). The axis of the H2 jet has a position angle of ~ 



23°. The H2 jet is almost in the plane of the sky, with the 
blueshifted side to the north and the redshifted side to the 
south of the source. 

3.1. 1.33 mm Continuum Emission 

Continuum emission is detected around the source at 
A = 1.33 mm with a flux of 0.11±0.03 Jy in between the 
diffuse nebulae seen in H2 (Fig. 1). This flux is simi- 
lar to the peak flux in an 11" beam around the source 
found by Chini et al. (1997) with the IRAM 30 telescope. 
The emission is not resolved and shows a slight elonga- 
tion along the jet axis. Previously, continuum emission 
has been detected around the source at A = 1.1 mm by 
the JCMT with a flux of 0.19 Jy (Zinnecker et al. 1992). 
Continuum emission has also been detected in the infrared 
by the IRAS. The fluxes corrected for extinction are 0.37, 
19.9, and 61.9 Jy at A = 25, 60, and 100 /xm, respectively. 
It was not detected at A = 12 jim. The spectral energy 
distribution (SED, see Fig. lb) is consistent with a Class 
source. It clearly suggests that the continuum emission at 
A = 1.33 mm is thermal emission from a dust core around 
the source. 

The temperature of the dust core can be derived from 
the SED. For simplicity, we assume a constant tempera- 
ture for the dust core, and the solid angle subtended by 
the dust core, fi, is independent of frequency, v, so that 
the observed flux is given by 



F^ = n{l-e-^'')B,{Taust) 
where the optical depth 
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with Kn being the mass opac;ity and a being the column 
density. With (see, e.g., Beckwith et al. 1990) 
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the optical depth can be simplified as 

T. = (4) 

with an optical depth of unity at Vq- As can be seen from 
the SED, the emission are optically thin at A = 1.1 and 
1.33 mm, and their flux ratio implies that (3^1. The 
solid angle Vt is uncertain, likely < 2.3" x 1" (see Fig. la). 
With Q. = 1.2 and 2.3 arcsec^, the SED can be fitted with 
Tdust = 47.6 K and = 3 X lO^^ Hz, and Tdust = 45.4 K 
and Vo = Q y- 10^2 Hz, respectively (see Fig. lb). Thus, 
Tdusi is assumed to be 46 K. With this temperature and 
the optically thin emission at A = 1.33 mm, the (gas -|- 
dust) mass is estimated to be 



M 



O.OSMp 







where the distance to the source D = 460 pc. 
3.2. 



(5) 



C^^O (J = 2-1) emission 

In C^^O, an envelope is seen around the source in the 
equator between the diffuse nebulae (Fig. 2a). Faint shells 

(2 (T contours) are also seen opening to the north and south 
surrounding the bases of the nebulae. The envelope is bet- 
ter seen in the blueshifted emission, in which the shells are 
much fainter than the envelope (Fig. 2b). It is extended 
across the source from the east to the west with a position 
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angle of ^ 110°±10°, perpendicular to the jet axis. It is 
flattened with an (deconvolved) FWHM extent of ~ 4.3" 
(or 2000 AU at 460 pc) and an aspect ratio of ~ 3:2. In the 
redshifted emission, the envelope is less extended. Since 
the redshifted emission around the bases of the nebulae is 
as bright as that of the envelope, the redshifted emission 
appears extended along the jet axis. 

A cut along the major axis of the flattened envelope 
shows a PV structure with the emission seen across the 
soiirce with the redshifted and blucshiftcd peaks on ei- 
ther side of the source (Fig. 3). The blueshifted emission 
spreads across the source with a peak at (—0.5 km s~^, 
—0.45") in the west, while the redshifted emission spreads 
across the source with a peak at (0.35 km s~^, 0.1") in 
the east. This PV structure is similar to that seen in 
C^^O J=l-0 toward an edge-on dynamically infalling en- 
velope with rotation around a young Class source IRAS 
04368+2557 in L1527 (Ohashi ct al. 1997). It suggests 
that the material of the flattened envelope is also infalling 
toward the source with rotation, with the redshifted emis- 
sion from the nearside, while the blueshifted emission from 
the farside. 

3.3. ^^CO (J = 2-1) emission 

In ^^CO, the emission peaks at the source and extends 
to the bases of the nebulae (Fig. 4a). The emission can 
be decomposed into low-velocity and high-velocity com- 
ponents (sec Fig. 5b and the explanation in the following 
paragraph). At low velocity, the emission forms a hour- 
glass structure with a narrow waist around the source (Fig. 
4b), with the blueshifted emission opening to the north 
and the redshifted emission opening to the south around 
the bases of the nebulae (Fig. 4c). Around the source, 
the blueshifted emission peaks a little to the west, while 
the redshifted emission shifts to the east of the source and 
is weak along the major axis of the flattened C^^O enve- 
lope. At high velocity, the blueshifted emission peaks to 
the west while the redshifted emission peaks to the east 
of the source (Fig. 4d). The emission also extends to the 
north and south along the jet axis. 

A cut across the ^''CO waist along the major axis of the 
flattened C^^O envelope shows that the ^^CO emission 
has a larger velocity range than that of the C^^O emission 
(Figs. 5a and 5b). The velocity range is thus decomposed 
into low-velocity and high-velocity components, with the 
range of the low- velocity component similar to that of the 
C"^^0 emission. At low velocity, the blueshifted emission 
is seen across the source along the major axis of the flat- 
tened envelope with a peak at (—0.3 km s~^, 0.1"), and a 
redshifted dip is seen around 0.4 km s~^ across the source 
along the major axis of the flattened envelope. This dip 
may be due to the absorption of a cold material infalling 
toward the source (Evans 1999), which is probably traced 
by the bright C^^O emission around that velocity (see Fig. 
5b) . At high velocity, the redshifted and blueshifted emis- 
sion are each confined on one side of the source, as seen 
in Figure 4d. This high-velocity emission is unlikely to 
be from outflow, with the redshifted and blueshifted emis- 
sion on either side of the source. This high-velocity emis- 
sion may trace the inner part of the envelope, where the 
rotation starts to dominate over the infall motion and a 
rotationally supported disk starts to form. 



3.4. CO (J = 2-1) emission 

A ^^CO outflow is seen surrounding the H2 jet, with 
a narrow waist around the source (Fig. 6). Near the 
source, limb-brightened outflow shells are seen opening to 
the north and south from the source surrounding the H2 
jet. The shells may extend further down along the jet axis, 
connecting to the prominent II2 bow shocks SB4 and NB3, 
as seen in lower J transition line of ^^CO (Lee et al. 2000). 
Deeper observations with shorter baselines are needed to 
confirm this. At high velocities, the ^^CO emission is jet- 
like, i.e., coUimatcd and knotty, much like the H2 jet. lo- 
cated inside the cavity defined by the diffuse nebulae and 
i^CO outflow shells (Fig. 6b). In the south, ^^CO outflow 
shell is seen coincident with the H2 bow shock SB4, proba- 
bly tracing the ambient material swept up by the H2 bow 
shock. Note that here the ^^CO outflow shell is mainly 
seen on one side of the bow shock, which is brighter in H2. 

Assuming an optically thin emission, a temperature of 
50 K, and a ^^CO abundance of 8.5xl0~^ with respect 
to molecular hydrogen (Frerking et al. 1982), the H2 col- 
umn density of the shells around the source and the bow 
shock SB4 is found to be ~ 10^° cm~^. Assuming that 
the shell thickness is ~ 3" (see Fig. 6), the density is ~ 
5x10'^ cm^"^. The peak column density along the jet axis 
is ~ 5 X 10^° cm^^. Assuming that the thickness is ~ 
2", the density is ~ 4x10'' cm"^. Note that due to the 
absorption of the ambient cloud, that part of the emission 
is optically thick, and that part of the emission is resolved 
out by the interferometer, the values presented here are 
lower limits of the true values. 

This outflow, which is almost in the plane of the sky, is 
one of the best candidates to study the outflow kinematics 
transverse to the jet axis. A cut along the jet axis shows 
a series of convex spur PV structures on both redshifted 
and blueshifted sides with the highest velocity near the H2 
bow tips and knots (Fig. 7). Notice that the redshifted 
emission around ~ 7 km s~^ merges with the foreground 
ambient cloud emission in velocity and thus is mostly re- 
solved out from our observations. These PV structures 
are similar to those seen in lower J transition line of ^^CO 
(Lee et al. 2000), but with a higher velocity around the 
H2 knots and bow tips. These PV structures suggest that 
the outflow is accelerated by the shocks localized at the 
H2 bow tips and knots, so that the velocity of the outflow 
decreases rapidly away from the tips, similar to that seen 
in a pulsed jet simulation (see, e.g., Lee et al. 2001). 

3.5. SO (Jk = 65 — ^i) emission 

SO emission is detected around the source and the H2 
knots in the south, especially the knot SK4 (Fig. 8a). 
Around the source, the emission is not well resolved and 
is slightly elongated with a major axis about mid-way be- 
tween the jet and envelope axes. The blueshifted emis- 
sion is to the north, while the redshifted emission is to 
the south of the source (Fig. 8b), similar to that seen 
in the H2 jet, suggesting that the SO emission may arise 
from interactions with the jet. The peaks of the redshifted 
and blueshifted emission, however, are a little off the jet 
axis, with the redshifted one shifted to the east and the 
blueshifted one shifted to the west, suggesting that the 
SO emission may also arise from the envelope. Around 
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the knot SK4, the redshifted emission is to the east while 
the blueshifted emission is to the west, similar to that seen 
in the flattened C^^O envelope. 

The SO emission around the source may consist of more 
than one component. A cut along the jet axis shows two 
peaks at high velocity, one is blueshifted at (—5.5 km s~^, 
0.64") in the north and one is redshifted at (6.1 km s~^, 
—0.53") in the south (Fig. 9a). The rapid increase in the 
velocity suggests that the SO emission at high velocity in- 
deed arisc's from shock interactions with the jet. However, 
the SO emission seen around the systemic velocity may not 
be produced by the shock interactions, because of its low 
velocity. A cut across the jet axis shows a rather symmet- 
ric PV structure about the source (Fig. 9b) at low veloc- 
ity. The spatial extent and the velocity distribution (i.e., 
the redshifted emission is to the east and the blueshifted 
emission is to the west) of this PV structure are similar to 
that of the flattened C-'^^O envelope, suggesting that the 
SO emission at low velocity may be associated with the 
flattened envelope. There are also high blueshifted and 
redshifted emission at ^ 3" in the east. 

A cut across SK4 shows that the redshifted emission is 
to the cast, while the blueshifted emission is to the west 
(Fig. 9c), as seen in Figure 8b. The velocity extent seems 
to be similar to that around the source, with the higher 
velocity near the jet axis. 

3.6. Summary 

In this section, we summarize our observations with two 
composite figures, Figures 10a and 10b. From these fig- 
ures, it is clear that (1) the dust core seen in continuum 
and the -'^'^CO waist are surrounded by the flattened C^^O 
envelope in the equator, and thus may both trace the in- 
ner part of the envelope; (2) the ^''CO emission extends 
to the north and south from the source, surrounding the 
bases of the nebulae; (3) the ^^CO outflow appears as limb- 
brightened shells surrounding the H2 jet, and (4) the SO 
emission is detected around the source and the H2 knots 
in the south. 

4. MODELING OF CIRCUMSTELLAR ENVELOPE 

In the following, a simple edge-on disk model with both 

infall and rotation is used to reproduce the structures and 
kinematics of the flattened C^*0 envelope and the ^'^CO 
waist seen around the source. In this model, the disk is 
assimied to have a constant thickness of H, an inner ra- 
dius of i?in, and an outer radius of Rout- The number 
density of molecular hydrogen is assumed to be given by 
(in Cylindrical coordinates) 

n = no{^r (6) 

where no is the density at a characteristic radius of Rq = 
460 AU (i.e., 1") and p is a power- law index. The tem- 
perature of the envelope is uncertain and assumed to be 
given by 

T = To{^y (7) 

where Tq is the temperature at i?o and g is a power-law in- 
dex. The abundances of C^^O and ^'^CO relative to molec- 
ular hydrogen are assumed to be constant and given by 



1.7 X 10 ^ and 9.5 x 10 ^, respectively (Frerking et al. 
1982). 

The detailed radial profiles of the infall velocity and the 
rotation velocity can not be determined from our observa- 
tions obtained at current spatial and velocity resolutions. 
For simplicity, we assume a dynamical infall that results in 
a free fall for the infall velocity. Assuming that the mass 
of the envelope is small compared with that of the source, 
we have 

/ 2GM7 / M. Rq , _i 

"«=V^^'-''V0T5Aio^ 
where is the mass of the source. The rotation seems 
differential with the velocity increasing toward the source 
(sec Figs. 3 and 5). However, the rotation is unlikely to 
be Keplerian, because the rotation velocity, as found in 
the following, is much smaller than the infall velocity. In 
a dynamically infalling envelope with slow rotation, spe- 
cific angular momentum of each gas element is considered 
to be conserved, until the infall motion shifts to the ccn- 
trifugally supported motion around the radius where the 
rotation velocity is comparable to the infall velocity (see, 
e.g., Nakamura 2000). Thus, the rotation velocity is as- 
sumed to be given by 

V0 = V4,o— (9) 

where v^o is the rotation velocity at Rq and it depends on 
the initial angular momentum at the outer radius. In the 
model calculations, radiative transfer is used to calculate 
the emission, with an assumption of local thermal equilib- 
rium. For simplicity, the line width is assumed to be given 
by the thermal line width only. The line width due to tur- 
bulence is not included in our model. The channel maps 
of the emission derived from the model are convolved with 
the observed beams and velocity resolutions, and then used 
to make the integrated maps and PV diagrams. 

There are nine parameters in this model: H, i?in, i?out, 
r?,o, p. To, q, M,, and v^q. To narrow our search, only 
two promising values oi p are considered: (1) p = —1.5, 
as found in many theoretical infalling models (see, e.g., 
Nakamura 2000) and HCO"*" observations around low- mass 
envelopes (Hogcrhcijde 2001), and (2) p — —2.0, as found 
in 850 and 450 /im continuum observations of low-mass 
envelopes (Shirley et al. 2000). We find that these values 
of p can both result in a reasonable fit to our observations 
too (sec Fig. 11). In addition, the values of other param- 
eters do not depend much on these values of p. In our 
fits, we have H - 550 AU (- 1.2"), i?in ~ 115 AU (~ 
0.25"), Ront ~ 2300 AU (~ 5"), Af, - O.15M0, v^o ~ 0.3 
km , and no ~ 4 x 10^ cm~^. To and g, however, are 
required to be different between C^*0 and ^"'CO in order 
to fit the observations. We have To 35 K and q ~ —0.4 

for C^^O, and Tq ~ 23 K and q 1 for ^^CO. With 

these values of parameters, our model can roughly repro- 
duce: (1) the observed structures of the flattened C^^O 
envelope and the ^'^CO waist seen around the source; (2) 
the observed PV structures with the emission seen across 
the source with the redshifted and blueshifted peaks on 
either side of the source; and (3) the observed redshifted 
dip in ^■'CO. Note that in the observations, the emission 
extending to the north and south is likely from the enve- 
lope material affected by the outflow shells (see §5.1) and 
thus is not included in our model. 
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The mass of the flattened envelope, which is given by 

Me = 2mnnaH / {—f2nRdR (10) 

is ~ 0.028 and 0.025 M© for p = -1.5 and -2, respec- 
tively. This mass, however, is a factor of 3 lower than that 
of the dust core. Note that the mass derived for the dust 
core could have an uncertainty with a factor of 5 due to 
the large uncertainty in the mass opacity (Beckwith et al. 
1990). The infall rate, which is given by 

M{R) = 2'KRHn\vR\2mB. (11) 
is ~ 6 X 1O-*^M0 yr-i at Ro = 460 AU. 

5. DISCUSSION 

HH 212 is a highly coUimatcd jet powered by a young 
Class source, IRAS 05413-0104. We have mapped the jet 
and its source in 1.33 mm continuum, -"^^CO, ^"^CO, C^*0, 
and SO. In the following, we discuss our results in details. 

5.1. Circumstellar Envelope 

A flattened envelope is seen in C^^O around the source 
in the equator perpendicular to the jet axis, with its inner 
part seen in ^^CO. The structure and kinematics of the en- 
velope can be roughly reproduced by a simple edge-on disk 
model with both infall and rotation. In this model, the 
density of the disk is assumed to have a power-law index 
of p = —1.5 or —2, as found in other low- mass envelopes 
(Shirley ct al. 2000; Hogerheijdc 2001). The temperature 
of C^^O has a power-law index of g ~ —0.4, as found in 
other low-mass envelopes (Shirley et al. 2000; Hogerheijde 
2001). However, in order to reproduce the compact waist 
and the redshifted dip in ^"^CO, the temperature of ^'^CO 
is required to decrease with g ~ — 1, which is faster than 
that of C^^O. Notice that in our model, the redshifted dip 
is assumed to be produced by the absorption of a cold 
material at the outer part of the envelope. However, the 
dip may be partly due to missing flux in our interferomet- 
ric observations (Gueth et al. 1997). Further observations 
and modeling (e.g., with flared disk model, turbulence, and 
missing flux consideration) are needed to improve this. 

The envelope seems dynamically infalling toward the 
source with slow rotation because the kinematics is found 
to be roughly consistent with a free fall toward the source 
plus a rotation of a constant specific angular momentum. 
If this is the case, the source has a mass of ^ 0.15 Mq 
and the infall rate is ~ 6 x IO'^Mq yr^^ at Rq = 460 AU, 
similar to those found in other flattened infalling envelopes 
with slow rotation around low-mass protostars (Ohashi et 
al. 1997; Momosc et al. 1998). Assuming that the accretion 
rate is constant in the past and given by the infall rate, the 
accretion time will be M*/M ~ 3 x 10^ yr, as expected 
for a Class source. Notice that in theoretical models 
(e.g., Shu et al. 2000), this accretion rate corresponds to 
an isothermal sound speed of 0.3 km s""'^, or equivalently 
a temperature of 25 K, which is consistent with our simple 
model. This dynamically infalling envelope with rotation 
is expected to form a rotationally supported disk at i? ^ 74 
AU 0.16"), at which the rotation velocity is compara- 
ble to the infall velocity (Hayashi et al. 1993; Lin et al. 
1994). The ^^CO emission seen at high velocity may arise 
from this rotationally supported disk, with the redshifted 



emission and the blueshifted emission on cither side of the 
source (Fig. 4d). 

A flattened NH3 envelope has been seen rotating around 
the source in the same direction (Wiseman et al. 2001). It 
has a characteristic radius of ^ 3500 AU (^ 8") and can be 
considered as the extension of the flattened C^^O envelope 
into the surrounding. At that radius, the infall velocity 
and rotation velocity are expected to be ~ 0.27 and 0.04 
km s~^, respectively in our simple model. This rotation 
velocity is consistent witli that observed in the NH3 en- 
velope, suggesting that the rotation of the C^^O envelope 
is connected to that of the NH3 envelope, and that the 
angular momentum is carried inward into the C^^O enve- 
lope from the NH3 envelope. However, no infall velocity 
has been seen in the NH3 envelope. In theoretical models 
(e.g., Shu 1977), the infall radius is given by the isothermal 
sound speed times the accretion time. Assuming a tem- 
perature of 25 K, the isothermal soimd speed is ^ 0.3 km 
s""'^, resulting in an infall radius of ~ 2000 AU, similar to 
the outer radius of the C^^O envelope. Thus, it is possible 
that no infall motion is seen in the NH3 envelope. 

The flattened envelope is not formed by rotation be- 
cause it is not rotationally supported, with the rotation 
velocity smaller than the infall velocity. Could it result 
from an interaction with the CO outflow? The flattened 
NH3 envelope is seen carved by the CO outflow, forming 
a bowl structure around it (Wiseman et al. 2001; Lee et 
al. 2000). In our observations, the shells seen in C^*^0 
and ^^CO around the bases of the nebulae may trace the 
(swept-up) envelope material pushed by the outflow shells 
toward the equator, suggesting that the envelope around 
the source is also carved by the outflow. Having said that, 
the flattened envelope may not form by the outflow in- 
teraction, because the shells are seen separated from the 
envelope. 

Formation mechanisms of such a flattened infalling en- 
velope around a forming star have been proposed with 
and without magnetic field. Galli & Shu (1993) examined 
the gravitational contraction of magnetized spherical cloud 
cores. They showed that flattened infalling envelopes are 
formed as the magnetic field impedes the contraction per- 
pendicular to the field lines, which are aligned with the 
rotation axis. Hartmann et al. (1996) showed that even 
in the absence of a magnetic field, initially sheetlike cloud 
cores can also collapse to form flattened envelopes. Re- 
cently, Nakamura (2000) also showed that flattened in- 
falling envelopes are a natural outcome of the gravitational 
contraction of prolate cloud cores with slow rotation. In 
this scenario, the flattened envelopes are predicted to have 
n oc r^^ "'', Vr oc r^'^'^, and oc r^^, the same as our 
model. 

5.2. Molecular Outflow 

A ^^CO outflow is seen surrounding the H2 jet, with a 
narrow waist around the source. The outflow seen here in 
12C0 J = 2 - 1 is similar to that seen in ^^CO J = 1 - 0, 
but with higher J transition line of ^^CO tracing higher ve- 
locity around the knots and bow tips. The morphological 
relationship between the H2 jet and the ^^CO outflow, and 
the kinematics of the ^^CO outflow along the jet axis are 
both consistent with those seen in a jet-driven bow shock 
model (see, e.g., Lee et al. 2001, and reference therein). 
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In that model, a highly coUimatcd jet is launched from 
around the source, propagating into the ambient medium. 
A temporal variation in the jet velocity produces a chain 
of H2 knots propagating down along the jet axis. Since 
these H2 knots are localized (ballistic) shocks with high 
thermal pressure, they expands sideways, growing into H2 
bow shocks with the shock velocity (and thus tempera- 
ture) decreasing rapidly away from the bow tips. As these 
bow shocks propagate down along the jet axis, they sweep 
up tlic^ ambient material into a thin outflow shell around 
the jet axis. In our observations, the ^^CO shells around 
the source can be identified as the ambient material swept 
up by the II2 bow shocks. In addition, ^^CO emission also 
arises from aroimd the II2 knots and bow shocks, with the 
velocity decreasing rapidly away from the knots and bow 
tips. Higher J transition line of ^^CO traces higher veloc- 
ity around the knots and bow tips likely because it arises 
from regions of higher temperature and thus higher veloc- 
ity closer to the knots and bow tips. As a result, the ^^CO 
emission appears jetlike at high velocities in high J tran- 
sition. Located inside the outflow cavity, the jetlike ^^CO 
emission is unlikely from the entrained ambient gas. The 
jet is probably intrinsically molecular, such as that seen 
in HH 211 (Gueth & Guilioteau 1999). Since the shell is 
the swept-up ambient material, the density of the ambi- 
ent material is expected to be lower than that of the shell. 
Thus, the jet is an overdensed jet with the density higher 
than that of the ambient material. 

In the south, ^^CO emission is seen coincident with the 
H2 bow shock SB4, indicating that the shocked material 
there must have cooled very fast from ^ 2000 K (H2) to a 
few 10 K (^^CO) because of a significant mixing with the 
ambient material. The transverse velocity there is also 
lower than that around other bow shocks (see Fig. 7), 
probably because the bow shock SB4 is sharing momen- 
tum with the quiescent ambient material (cf. Fig. 7 in 
Lee et al. 2001). Therefore, the bow shock SB4 may trace 
the leading bow shock where the head of the jet impacts 
on the ambient material. 

5.3. SO emission 

SO emission is seen around the source and the II2 knotty 
shocks in the south, especially the knot SK4. The emis- 
sion around the source may arise from two components. 
Like that around the knot SK4, the emission seen at high 
velocity is likely to be from knotty shocks, which are ex- 
pected to be there following the H2 knotty shocks toward 
the source. The shocks are not seen in H2 probably due 
to dust extinction. However, the SO emission seen around 
the systemic velocity may arise from the flattened enve- 
lope due to its low velocity. Observations at high angular 
resolution are needed to confirm this. 

SO emission has also been seen toward the CepA- 
East outflows, tracing ambient material at 60-100 K and 
shocked gas at 70-180 K (Codclla et al. 2005). Assuming 
an optically thin emission, a temperature of 120 K, and 
a SO abundance of 10~^ relative to molecular hydrogen 
(Codclla et al. 2005), the mass is found to be ~ 0.4 and 
0.1 Mq, respectively, around the source and the knot SK4. 
However, these masses arc too large to be from knotty 
shocks. The abundance of SO must have been highly en- 
hanced. Note that ^^CO emission is weak around the knot 



SK4 (see Fig. 6b), probably suggesting that ^^CO there is 
destroyed by the shock interaction. The detection of the 
SO emission is likely because of the abundance enhance- 
ment of SO due to the shock-triggered release of various 
molecules (e.g., H2S/OCS/S) from dust mantles and the 
subsequent chemical reactions in the warm gas. Since the 
abundance is enhanced only for a very short period of time, 
the SO emission traces regions of recent shock activity 
(Codella et al. 2005). 

Around the knot SK4, the redshiftcd emission is seen to 
the east while the blueshifted emission is seen to the west, 
similar to that seen in the flattened C^^O envelope. A sim- 
ilar velocity structure has been seen around the knot SKI 
(sec Fig. 8 for its location) in H2 and considered as a ten- 
tative evidence that the knot SKI is rotating (Davis et al. 
2000). It was observed with 3 slits from east to west with 
a separation of ~ 0.45". The peak velocities with respect 
to the systemic velocity were found to be +4.6, +2.9, and 
+2.3 km s^^, measured from east to west. Thus, assum- 
ing that the central velocity at the knot SKI is 3 km 
s~^, the redshifted emission is seen to the east while the 
blueshifted is seen to the west. A similar velocity struc- 
ture is seen here around the knot SK4 in SO but with 
higher velocity. Does this suggest that the knot SK4 is 
also rotating and that the jet carries away angular mo- 
mentum from the infalling envelope? The velocity here, 
however, is much larger than v^q (i.e., 0.3 km s~^) and 
thus may be too large to be from rotation. In addition, 
the shocked material there is expected to have sideways 
expansion. Further observations at higher angular resolu- 
tion are needed to study it. 

6. CONCLUSION 

We have mapped the 1.33 mm continuum, ^^CO (J = 
2 - 1), 13C0 (J = 2 - 1), C^^O (J = 2 - 1), and SO 
{Jk = 65 — 54) emission around a protostellar jet HH 212 
and its central source, IRAS 05413-0104. A dust core is 
seen in the continuum around the source, with a temper- 
ature of ^ 46 K and a mass of ^ 0.08 Mq. A flattened 
envelope is seen in C^^O around the source in the equator 
perpendicular to the jet axis, with its inner part seen in 
^'^CO. The structure and kinematics of the envelope can be 
roughly reproduced by a simple edge-on disk model with 
both infall and rotation. The flattened envelope seems dy- 
namically infalling toward the source with slow rotation 
because the kinematics is found to be roughly consistent 
with a free fall toward the source plus a rotation of a con- 
stant specific angular momentum. A ^^CO outflow is seen 
surrounding the H2 jet, with a narrow waist around the 
source. Jctlike structures are also seen in ^^CO near the 
source aligned with the H2 jet at high velocities. The mor- 
phological relationship between the H2 jet and the ^^CO 
outflow, and the kinematics of the ^^CO outflow along the 
jet axis are both consistent with those seen in a jet-driven 
bow shock model. SO emission is seen around the source 
and the H2 knotty shocks in the south, tracing shocked 
emission around them. 
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and maintaining the array. H.B. acknowledges financial 
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Forschungsgemeinschaft (DFG, grant BE2578/1). 
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Table 1 
Summary of Observations 



Line 


Frequency 


Beam Size 


Channel Width 


Noise Level" a 




(GHz) 


(arcsec x arcsec) 


(km s""'^) 


(Jy beam"-'^) 


1.33 mm continuum 


225.000000 


2.0x2.3 




3.5E-3 


12C0 J = 2 - 1 


230.537980 


2.8x2.3 


0.264 


0.25 


13C0 J = 2 - 1 


220.398676 


2.8x2.4 


0.274 


0.15 


C^^O J = 2 - 1 


219.560357 


2.8x2.4 


0.274 


0.14 


SO Jk = 65 — 54 


219.949442 


2.8x2.4 


0.554 


0.12 



''RMS noise level per velocity channel 



8 




Aa (") V (Hz) 

Fig. 1. — (a) 1.33 mm continuum contours plotted on top of the H2 image adopted from McCaughrean et al. (2002). The contours go 
from 3 to 18 O" with a step of 3 cr, where a = 3.5 mjy beam~^. The asterisk indicates the IRAS source. Note that the (0, 0) position here 
corresponds to a(2ooo) = 05''43™51=30, 5(2000) = -01°02'53.0". (6) Simple fits to the BED of the continuum source with Q = 1.2 (solid Une) 
and 2.3 (dashed line) arcsec^. 




10 5 -5 10 5 -5 

A« (") 

Fig. 2. — C^^O emission plotted on top of the H2 image, with the asterisk indicating the IRAS source, (a) Total emission integrated over 
2.2 km s~^ from —1.1 km s~^ to 1.1 km s~^. The contours go from 2 to 23 cr with a step of 3 a, where a = 0.11 Jy beam~^. (6) Redshifted 
(red contours) and blueshifted (blue contours) emission, integrated over 1.1 km s~^ on the redshifted and blueshifted side, respectively. The 
contours go from 2 to 17 cr with a step of 3 cr, where cr = 0.08 Jy beam~^ km s~^. 
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Fig. 3.— PV diagram of C^^O emission cut along the major axis of the flattened envelope. 
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Fig. 4. 



^CO emission plotted on top of the H2 image, with the asterisk indicating the IRAS source, (a) Total emission integrated over 



3.4 km s~^ from -1.7 km s~^ to 1.7 km s~^. The contours go from 2 to 20 cr with a step of 3 cr, where a = 0.15 Jy beam~^ km s~^. (b) 
Low-velocity emission integrated from -0.7 to 0.7 km s~^. The contours go from 2 to 17 cr with a step of 3 cr, where a = 0.09 Jy beam^^ km 
s~^. (c) Low-blueshifted (blue contours) and low-redshifted (red contours) emission, integrated from -0.7 to km s~^ and to 0.7 km s~^, 
respectively. The contours go from 2 to 18 (t with a step of 2 a, where a = 0.06 Jy beam~^ km s~^. (d) High-blueshifted (blue contours) 
and high-redshifted (red contours) emission, integrated from -1.7 to -0.7 km and 0.7 to 1.7 km s~^, respectively. The contours go from 2 
to 14 (T with a step of 2 cr, where cr = 0.08 Jy beam^^ km s^^. 



11 




Fig. 5. — PV diagrams of C^*^0 and ^"^CO emission cut along the major axis of the flattened envelope, (a) for ^"^CO. (b) Overlay of the 
two. Here, LV denotes low velocity and HV denotes high velocity. 
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Fig. 6. — -'^■^CO emission plotted on top of the H2 image adopted from McCaughrean et al. (2002), with the asterisk indicating the IRAS 
source. The black dotted line outlines the observed region in our observations, (a) Black contours: Total emission integrated over 32 km s^^ 
from —16 to 16 km s^^. The contours go from 2 to 24 ct with a step of 2 cr, where ct = 0.68 Jy bcam^-"^ km s^^. (6) Black contours: High 
blucshifted emission integrated from —5 to —16 km s^^. The contours go from 3 to 21 a with a step of 3 ct, where a = 0.42 ,Iy bcam^-'^ km 
s^'^. Gray contours: High redshifted emission integrated from 2 to 13 km s~^. The contours go from 3 to 12 cr with a step of 3 it, where 
u = 0.42 Jy beam^^ km s~^. Here, the velocity ranges are selected to show the jetlike structures as much as possible, without the confusion 
from the shell emission. 
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V , (km/ s) 
obs ^ ' ' 

Fig. 7. — Left: The H2 jet and CO outflow, with the outflow axis rotated to be aligned with the y-axis, The two vertical lines indicate the 
width of the cut for the position-velocity diagram. Right: Position-velocity diagram of the CO outflow cut along the outflow axis. The two 
dotted lines around 7 km s^^ indicate the velocity range of a foreground ambient cloud seen by the FCRAO single-dish telescope (Lee et al. 
2000). 
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Fig. 8. — SO emission plotted on top of the H2 image, with the asterislt indicating the IRAS source, (a) Total emission integrated over 16 
km from -7 to 9 km s~^. The contours go from 3 to 13 it with a step of 2 a, where a = 0.35 Jy beam~^ km s~^. (6) Redshifted (red 
contours) and blueshifted (blue contours) emission, integrated from -7 to -2 km s~^ and from 2 to 9 km , respectively. The contours go 
from 3 to 10 cr with a step of 1 a, where a = 0.23 Jy beam~^ km s~^. 




Fig. 9. — PV diagrams of SO emission, (a) A cut along the jet axis for the emission around the source. (6) A cut across the jet axis for 
the emission around the source. LV denotes low velocity, (c) A cut across the jet axis for the emission around the knot SK4. 
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Fig. 11. — Modeling of envelope in C-^^O {a,b,e,f) and -'^''CO {c,d,g,h). Gray eontours with image are the observations, blaek eontours 
are derived from our simple model. (a)-{d) shows the total emission, while (e)-(/i) shows the PV diagrams cut along the major axis of the 
flattened envelope. 



